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I. INTRODUCTION 



The flavor-changing neutral current transitions B — > K + X and B c — > D(D*) + X with X = 7, are of 

special interest because they proceed at the loop level in the Standard Model (SM) involving also the top quark. 
They may therefore be used for a determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vt q 
(q = d,s,b). The available experimental measurements of the branching ratio of the inclusive radiative B-meson 
decay 

' (3.11 ± 0.80(stat) ± 0.72(syst)) x 1(T 4 ALEPH [1] 

Br (B -> X sl ) = \ (3.36 ± 0.53(stat) ± 0.42(syst)±°-j°(th)) x 1CT 4 BELLE [2j 
_ (3.21 ± 0.43(stat) ± 0.27(syst)±°;^(th)) x 1CT 4 CLEO [3] 
are consistent with the next-to-leading order prediction of the standard model (see, e.g. Q and references therein): 

Bt(B -> A s7 ) SM = (3.35 ± 0.30) x 10~ 4 . (1) 
The decay B — > K (I = e, fj,) has been observed by the BELLE Collaboration [j| with a branching ratio of 

Br (B -» K l+r) = (0.75±° H ± 0.09) x 10" 6 . (2) 

The recent observation of the bottom-charm B c meson by the CDF Collaboration at Tevatron in Fermilab @ raises 
hopes that one may also explore the rare decays of the bottom-charm meson in the future. 
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The theoretical study of the exclusive rare decays proceeds in two steps. First, the effective Hamiltonian for 
such transitions is derived by calculating the leading and next-to-leading loop diagrams in the SM and by using the 
operator product expansion and renormalization group techniques. The modern status of this part of the calculation 
is described in the review [?| (and references therein). Second, one needs to evaluate the matrix elements of the 
effective Hamiltonian between hadronic states. This part of the calculation is model dependent since it involves 
nonperturbative QCD. There are many papers on this subject. The decay rates, dilepton invariant mass spectra and 
the forward-backforward asymmetry in the decays B — > K (I = e, /i, r) have been investigated in the SM and its 
supersymmetric extensions by using improved form factors from light-cone QCD sum rules |8( ■ An updated analysis 
of these decays has been done in Q| by including explicit 0(a s ) and Aqcd/w& corrections. The invariant dilepton 
mass spectrum and the Dalitz plot for the decay B — > Kl + l~ have been studied in |j| by using quark model form 
factors. The B — > K decay form factors were studied via QCD sum rules in [ljj and within the lattice-constrained 
dispersion quark model in . Various aspects of these decays were discussed in numerous papers byAliev et al. pj. 
The exclusive semileptonic rare decays B — ► Kl + l~ were analyzed in supersymmetric theories in |l3| . The angular 
distribution and CP asymmetries in the decays B — ► Kne + e~ were investigated in ^4|. The lepton polarization for 
the inclusive decay B — > X s l + l~ was discussed in and The rare decays of B c — > D(D*)l + l~ were studied 
in [Tt| by using the form factors evaluated in the light front and constituent quark models. 

In this paper we study the exclusive rare decays B — > K 11. We employ a relativistic quark model 0,^1 to calculate 
the decay form factors. This model is based on an effective Lagrangian which describes the coupling of hadrons H 
to their constituent quarks. The coupling strength is determined by the compositeness condition Zfj — [2(3, [2l| 
where Zh is the wave function renormalization constant of the hadron H. One starts with an effective Lagrangian 
written down in terms of quark and hadron fields. Then, by using Feynman rules, the S-matrix elements describing 
the hadronic interactions are given in terms of a set of quark diagrams. In particular, the compositeness condition 
enables one to avoid a double counting of hadronic degrees of freedom. The approach is self-consistent and universally 
applicable. All calculations of physical observables are straightforward. The model has only a small set of adjustable 
parameters given by the values of the constituent quark masses and the scale parameters that define the size of the 
distribution of the constituent quarks inside a given hadron. The values of the fit parameters are within the window 
of expectations. 

The shape of the vertex functions and the quark propagators can in principle be found from an analysis of the Bcthe- 
Salpeter and Dyson-Schwinger equations as was done e.g. in 22]. In this paper, however, we choose a phenomenological 
approach where the vertex functions are modelled by a Gaussian form, the size parameter of which is determined by 
a fit to the leptonic and radiative decays of the lowest lying charm and bottom mesons. For the quark propagators we 
use the local representation. In the present calculations we do not employ the so-called impulse approximation used 
previously |l9j |. The numerical results obtained with and without the impulse approximation are close to each other 
for light-to- light and heavy-to-heavy transitions but differ considerably from one another for heavy-to-light transitions 
as e.g. in the B — ► 7r transitions. 

We calculate the form factors of the transition B — > K and use them to evaluate differential decay rates and polar- 
ization observables. We give the q 2 -dependence of a set of observables with and without long-distance contributions 
which include the lower-lying charmonium states according to [23] - We extend our analysis to the exclusive rare 
decay B c — > D(D*)ll. We derive a four-fold angular decay distribution for the cascade B c — > D*{-^ Dir)ll process 
in the helicity frame including lepton mass effects following the method outlined in |24j . The four- fold angular decay 
distribution allows one to define a number of physical observables which are amenable to measurement. We compare 
our results with the ones of other studies. 

We should remark that our approach is developed mainly for the hadrons (mesons and baryons) which satisfy to 
the so-called "threshold inequality": the hadron mass should be smaller the sum of their constituents, i.e. the sum 
of the constituent quark masses. In this vein, our model was successfully developed for the study of light hadrons 
(e.g., pion, kaon, baryon octet, A-resonance), heavy-light hadrons (e.g., D, D s , B and -B. s -mesons, Aq Ep, Eq and 
fig-baryons) and double heavy hadrons (e.g, J/*, T and £> c -mesons, Sqq and CIqq baryons) To extend 

our approach to other hadrons we had to introduce extra model parameters or do some approximations, like, e.g., 
to introduce the cutoff parameter for external hadron momenta to guarantee the fulfilment of the mentioned above 
"threshold inequality". Therefore, at the present stage we can not apply our approach for the study of rare decays 
involving K* mesons. Probably, it will be a subject of our future investigations. 

The layout of the paper is as follows. In Sec. II we discuss the effective Hamiltonian. We use the analytical 
expressions for the Wilson coefficients from and the input parameters from ||. Sec. Ill is devoted to the description 
of the B(B C ) — > K(D, D*)ll decays in terms of helicity amplitudes. We derive the four-fold angular decay distribution 
for the cascade B c — > D*(-^ Dtt)11 process in the helicity system and define a number of physical observables from 
the ang ular decay distribution. Our analysis goes beyond the results on the four-fold decay distribution presented in 
|25l 126) in that we include lepton mass effects appropriated for the treatment of the channel with r-leptons in the final 
state. We also present results on the longitudinal polarization of the leptons in the U-CM frame. This analysis differs 
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from the analysis in |25|, |26| where the polarization of the leptons were calculated in the initial meson rest system. 
We also include lepton mass effects in the polarization calculation. In Sec. IV we briefly discuss our relativistic 
quark model and demonstrate the difference between the exact calculation and the approximate calculation using the 
impulse approximation taking as an example the B — tt form factor. We calculate matrix elements and form factors 
for the decay B — > Kll and compare their behavior with those calculated in [8j- In Sec. V we present our numerical 
results for branching ratios and asymmetry parameters. We plot the g 2 -dependence of the differential decay rate and 
the longitudinal polarization of the leptons with and without long distance contributions. 



II. EFFECTIVE HAMILTONIAN 



The starting point of the description of the rare exclusive decays is the effective Hamiltonian obtained from the SM- 
diagrams by using the operator product expansion and renormalization group techniques. It allows one to separate 
the short-distance contributions and isolate them in the Wilson coefficients which can be studied systematically within 
perturbative QCD. The long-distance contributions are contained in the matrix elements of local operators. Contrary 
to the short-distance contributions the calculation of such matrix elements requires nonperturbative methods and is 
therefore model dependent. 

We will follow Refs. 7] in writing down the analytical expressions for the effective Hamiltonian and paper in 
using the numerical values of the input parameters characterizing the short-distance contributions. At the quark level, 
the rare semileptonic decay b — > s(d)l + l~ can be described in terms of the effective Hamiltonian: 

r io 

ff eff = --^A t ^C i ( M )Q i ( M ). (3) 

where A t = V^ d ^V t b is the product of CKM elements. For example, the standard set Q of local operators for 
b — > sl + l~ transition is written as 



(4) 



where G M „ and F^ v are the gluon and photon field strengths, respectively; Tjj are the generators of the SU(3) color 
group; i and j denote color indices (they are omitted in the color-singlet currents). Labels (V±A) stand for 7 M (1±7 5 ). 
Qi t 2 are current-current operators, Q3-6 are QCD penguin operators, Qt^s are "magnetic penguin" operators, and 
Qg^o are semileptonic electroweak penguin operators. Explicit formulae for the Wilson coefficients Ci(fi) obtained in 
leading logarithmic order are written down in Appendix A. 

The effective Hamiltonian leads to the free quark b — > sl + l _ -decay amplitude: 



Qi = 


(SiCj)v-A, (Cjbi)v-A, 


Q 2 = 


(sc) V - A (cb) V -A, 


Q3 = 


{sb) v ^ A Y,q{m)v-A, 


Qi = 


(sibj)v-A i2 q %q.i)v-A, 


Q 5 = 


(sb) V - A J2q(<l<l)V+A, 


Q 6 = 


(Sibj)v-A Y.q(qj1i)v+A, 


Qv = 


^m b (sa^(l+ l5 )b)F^, 


Qs = 


5 ^m 6 (5 i ff' u '(l + -ysjTyftjOG, 


Q9 = 


-^(sb)v-A(H)v, 


Qw - 


= -§^{sb) V -A{U)A 



2m b ^eff 



M(b -> sl + r) = Xt {Cf (sb) v ^ A (ll) v + Cio (sb) v „ A (II) A (5) 



2V^7T 

2 Cf (sia^(l + 7 5 )g i/ 6) (II), 



where Cf = C 7 - C 5 /3 - C 6 . 

The Wilson coefficient Cg ff effectively takes into account, first, the contributions from the four-quark operators Qi 
(i = 1, ...,6) and, second, the nonperturbative effects coming from the cc-resonance contributions which are as usual 
parametrized by a Breit-Wigner ansatz p3| : 

,..,.n ,-. \ 3-7T \ - TlVi — »■ l + l~) mvi 

Cf = C 9 + C < h(m c , s) + — k 2^ 



Vi=i/j(ls),ip(2s) 



i/i(M)(4C 3 + 4C4 + 3C 5 + C 6 ) (6) 
^/i(0, s) (C 3 + 3C 4 ) + I (3C 3 + C 4 + 3C 5 + C 6 ) . 
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where C* = 3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 . Here 

h(rh Cl s) = 



8, m h 8, . 8 4 

— In In m r H 1 — x 

9 u 9 27 9 



^(2 + z)|l-xr 



In 



2 arctan ■ 



r irni . -, 

for a; = — - < 1 



for x = > 1, 



h(Q,s) = 



27 



8 m b 4 4. 
— — in —in s + —nr. 



where m c 



mc/iriB, s — q 2 /m 2 B and n = 1/Cq. 



The relation between the MS b-quark mass rrib = m.{,(/it) at a scale y, and its pole mass mf, jPO i c is given by 



mb(/x) = TOfc^polo \ 1 - 



4 a s (p) 



l--ln 

4 



^&,polc 



(7) 



where 



-in 



1 - 



lnln-r# 
A 2 



QCD 



A) In 





1V QCD 



(3 2 o 



(8) 



hi 



A 2 

1V QCD 



and where /?o = 23/3 and /3i = 116/3, as is appropriate for five flavors. We will use a scale /i = mfc iPO i throughout 
this paper. The numerical values of the input parameters are taken from |8( and the corresponding values of the 
Wilson coefficients used in the numerical calculations are listed in Table [I] 



III. B -> Kll AND B c -» D(D*)?;-DECAYS 

A. Form factors and differential decay distributions 

We specify our choice of the momenta as p\ — p 2 + fci + k 2 with p\ = m\ , p\ = m 2 and k\ = k\ = y 2 where k\ 
and k 2 are the l + and Z~ momenta, and mi, m 2 , /i are the masses of initial meson H in , final meson ify and lepton, 
respectively. 

The matrix elements of the exclusive transitions B — > Kll and -B c — » D(D*)ll are defined by 

M(if in -> fl>Zl) = %-^-{C| ff <#/|sO"&| ff in > Z^i (9) 
V2 2 7r 

+ C 10 < HflsO^blHin > T-f-fl 
2rrib 



,2 



Cf < H f | at*"" (1 + 7 5 ) q v b\H m > Z 7 "Z 



where H m = B or B c , H f = K, D or D*. 
Wc define dimensionless form factors by 



< K(D)(p 2 )\s(d) lfl b\B(B c )( Pl ) >=F + (q 2 )P^+F^(q 2 )q^ (10) 

< K(D)(p 2 ) | s(d) ia^q" b I fl(S c )(Pi) >= z P m ^(9 2 ) . 

mi + m 2 p 



i<D*(p 2 ,e 2 )\d0^b\B c ( Pl ) >= \ e\ v 

mi + m 2 

x{-.9^Pg M<f ) + P^Pu A+{q 2 ) + q^P v A_(q 2 ) + ie^ aSi P a q Si V{q 2 )} . 



i < D*(p 2l e 2 )\di(j^q v {l + lb )b\B c (p 1 ) >= 

= 4 V { 9au Pqa (q 2 ) - Prt Pu a + (q 2 ) - ie^P^ g(q 2 )} 
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where P = pi + p 2 , q = Pi - P2, Pjt = p » - Q^Pq/q 2 , 9^ v = 9^ ~ g^/V, and 4 is tne polarization four-vector 
of the D*. Since we want to compare our calculations with those in [8| and (l7|. it is useful to relate our form factors 
to those used in 0] and 0|- The relations read 

F + = /Ali =F Gong 

F- = -^S^(/+-/o) Ali -^ Gcng 

TTi _ fAli _ jjiGcng 

= Af 1 = Al cng 



q 2 



V 


= v Kl{ 


_ _ yGeng 


do 


rpAYi 

— 1 2 


- o Gong 


a + 


= (T 2 - 


„2 

1 y > Ali _ „Gcng 
t" 2 2 J 3j a + 
— 


9 


- if 11 


= 5 Gons 



The matrix element in Eq jnj is written as 



Gf ocXt 



M (B(B C ) - A-(Z»*)«) = ^| ~ W (M + ^ OvY&O} (H) 

where the quantities Tf are expressed through the form factors and the Wilson coefficients in the following manner: 
(a) B(B C ) -> K{D)ll-decay: 

= P" + J*0> if (i = 1,2), (12) 

= CfF + + CfF T 2mb 



mi + m 2 



= C fF--CfF T -^-^ 



mi + 777,2 q 



= C W F ± . 



(b) B c -> L>*ZZ-decay: 



Tf = rr4, (i = l,2), (13) 



= {-Pqg^ Af + P^P V Af + q»P u A ( l ] + ie^ afi P a qa V (l) \ 

mi + 7772 

T/ (l) ^cff T/ , -off 2m b (mi +m 2 ) 

q 

4 1} = Cf A + Cfa 2mb{m y m2 \ 
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A« = Cf A + + Cf W^i + m,) 

= Cf A_ + Cf (oo - o+) 2Wfc(Wl 2 +W2) ^ , 

g z q A 

v (2) = CioV, 4 2) = c 10 a , a£ 2) = c 10 a±. 

Let us first consider the polar angle decay distribution differential in the momentum transfer squared q 2 . The polar 
angle is defined by the angle between q=pi—p2 and ki rest frame) as shown in Fig.Q] One has 

d 2 T \p 2 \v 1^,^,2 G% /a|A t |\ 2 |pa|t; 



dq 2 d cos 9 (27r) 3 4m3 8 ; 



I V IMI 2 = -^L ( ^ (14) 



x ^if 1 T-tr[7^(^-M)7,(^2 + /x)] 



+ ■ tr[7 M 7 5 (ft - fi) 7„75 (ft + /x)] 
+ ' tr ^ (#1 - M)7»75 (ft +M)] 



G| ^a|A t |\ 2 |p 2 |z; 1 



(2tt) 3 V 2vr y 8m? 2 \ "* 
- \ L$ ■ (q 2 Hg + (q 2 - 4 M 2 ) Hg ) + L$ • (^ + fl^)} 



where |p»2 1 = A 1 / 2 (m?, m 2 , q 2 )/2mi is the momentum of the K (D*)-meson and w = yl - 4/i 2 /g 2 is the lepton velocity 
both given in the i?(B c )-rest frame. We have introduced lepton and hadron tensors as 

(15) 

B. Helicity amplitudes and two-fold distributions 

The Lorentz contractions in Eq. (|14(l can be evaluated in terms of helicity amplitudes as described in First, 
we define an orthonormal and complete helicity basis e M (m) with the three spin 1 components orthogonal to the 
momentum transfer q' 1 , i.e. e AI (m)q M = for m = ±, 0, and the spin (time)-component m = t with e^(t) = q^ j \fq 2 ■ 

The orthonormality and completeness properties read 

e\(m)e ii (n)= g mn (m,n = i,±,0), 

(16) 

e^(m)el(n)g mn = 

with g mn = diag ( + , — , — , — ) . We include the time component polarization vector e M (t) in the set because we 
want to discuss lepton mass effects in the following. 

Using the completeness property we rewrite the contraction of the lepton and hadron tensors in Eq. 114(1 according 

to 

L^H% = L%e»\m)e^(m')g mm ,e^\ny(n')g nn ,H% 

= L^ k \m,n)g mml g nn ,H l:i {m',n') (17) 
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where we have introduced the lepton and hadron tensors in the space of the helicity components 

L^(m,n) = e?(m)J v (n)L$, H l3 (m,n) = e^(m)e v (n)H% . (18) 

The point is that the two tensors can be evaluated in two different Lorentz systems. The lepton tensors (to, n) will 
be evaluated in the ZZ-CM system whereas the hadron tensors 7J u (m, n) will be evaluated in the B(B C ) rest system. 
In the B(B C ) rest frame one has 

p£ - (mi, 0, 0, 0), 

p% = (£ 2 , 0, 0, -|pa|), (19) 
<f = (q Q , 0, 0, |p 2 |), 

where E% = (m\ + m\ — q 2 )/2mi and qo = (ml — m| + q 2 )/2mi. In the _B(B c )-rest frame the polarization vectors of 
the effective current read 

e"(t) = -^=(qo, 0, 0, |p 2 |), 



e"(±) - -±=(0, Tl, ~i, 0), (20) 
e"(Q) - -^=(|Pa|» 0, 0, q ) . 



Using this basis one can express the components of the hadronic tensors through the invariant form factors defined 
in Eq. (23J. 

(a) B(B C ) -> K(D) transition: 

H ij (m,n) = (et"(m)I^) • (e^{n)T^ = H l (m)H^(n) (21) 
The helicity form factors H l (m) are given in terms of the invariant form factors. One has 

H\t) = -±^(p q r + + q 2 r), 

H i (±) = 0, (22) 

fli (0 ) = 2mi ^r + . 



(b) B c — > 13* transition: 



H*(m,n) = et"(m)c"(n)fl^ = et"(m) C "(n)2^ (V + ^ W 

= e^(m)e>)T* a 4>)ef( S ),5 rs T^ (23) 
= £ t^(m)4 a (r)T; • (et>)e 2 ^( S )T4) V s = W{m)H^{n). 

From angular momentum conservation one has r = m and s = n for to, n — ±, and r, s = for to, n = t. For further 
evaluation one needs to specify the helicity components e 2 (m) (to = ±, 0) of the polarization vector of the D*. They 
read 

1 

(24) 



4(±) = — (0, ±1, -i, 0), 



e£(0) = — (IP2I, 0,0, -£ 2 

TO 2 
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They satisfy the orthonormality and completeness properties: 

e 2 A '( r ) e 2 M (s) = Srs, 

^(r)e\ y (s)5 rs = -g^ v 



(25) 



P2 l j,P2u 



'"-2 

Finally one obtains the non-zero components of the hadron tensors 

1 TOi|p 2 ' 



Wit) 


= ^(t)e\ a (Q)Tl a = - 


H\±) 


= e^{±)e\ a {±)T; a = 


£P(0) 


= ^(o) £ ^(o)r; Q 




= 1 v< 




mi + m 2 2m 2 y/q 2 



(Pq(-A a + A + ) + q 2 A^, 

mi + m 2 m 2 \/q 2 

i (-Pg4±2mi|p 2 |^ i ), (26) 

mi + 1712 



-Pq (m{ + m 2 - q 2 ) A a +Am{ |p 2 | 2 A\ 

The lepton tensors (m, n) are evaluated in the T/-CM system ki + k 2 — 0. One has (see Fig. ^) 

<f = (a/?, 0, 0, 0), 

fc^ 1 = (^i, |ki| sintfcosx, |k x | sinflsinx , |ki| cos 6» ) , (27) 
fcj = (Ei, — |kx| sin cos x , — |ki| sin # sin x , — |kx|cos#), 



with Ei — y/q 2 /2 and |ki| = \J q 2 — 4/i 2 /2. The longitudinal and time component polarization vectors in the 11 rest 
frame can be read off from Eq. I|20|l and are given by e M (0) = (0, 0, 0, 1) and e(t) = (1, 0, 0, 0) whereas the transverse 
parts remain unchanged from Eq. (|20|l . 

The differential (q 2 , cos 9) distribution finally reads 

dT(H m ^H f ll) 3 , 1/dT" dT 2 a 2 N 

dg 2 d(cos0) " 8 (1 + C ° S 6) '2 \-dcJ + ^) (28) 

3.2,1 fdri 1 dTf 
+ - sm 2 0.-1 —4- + — -§- 
4 2 V dg 2 dg 2 





3 , 
— u ■ - cos 6 
4 


dT 12 
dq 2 






+ ? sin 2 fl. 
4 








2 dg 2 




+ - cos 2 ■ 




3 . 2 
- - sin 

4 




2 dq 2 


obtains 








dr(ir m -> fl>S) 


1 


dr 22 


dr" 

+ 


dq 2 


2 ^ dq 2 


dq 2 


dg 2 






df 22 






2 dq 2 


dg 2 


2 dg 2 



cos 



dr 22 
2 ^ ai u 

dq 2 

dr 22 i df 22 

dg 2 4 dg 2 

Integrating over cos 9 one obtains 

dr 22 

" dg 2 

df 22 1 df 22 



(29) 



dg 2 2 dg 2 ' 

where the partial helicity rates dT^/dq 2 and df % /dq 2 (X = U, L, P, S; i, j — 1, 2) are defined as 

dTxy = G\ /a|A t |\ 2 |p 2 |g 2 ^ %3 
dq 2 (2tt) 3 V 2tt J 12 m 2 x ' 

dT Xij 2 y 2 dT% 
dq 2 q 2 dq 2 



(30) 
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The relevant bilinear combinations of the helicity amplitudes are defined in Table [H] 

To check our calculation we give the corresponding expression for the differential decay rate used in and : 



dr (B -> Kl+l-) _ G 2 F \\ t \ 2 m\a 2 



ds ~ 3 .2%5 VqjH 

and the expressions of 4>h, oii, f3, 
complete agreement with the decay distributions given in § and 



2t , 

1 + — <t> H a 2 + 12t#, 



s 



(31) 



with s = q 2 /m\ and t — fi 2 /m1 and the expressions of <j)jj, oti, /?», S (i=l,2,3) and 6 are given in Appendix B. We find 

17]- 



C. The four- fold angle distribution in the cascade decay 

B c ^D*(->Dn)Tl. 

The lepton-hadron correlation function L^H^ reveals even more structure when one uses the cascade decay 
B c — ► £)*(—> Dir)ll to analyze the polarization of the D*. The hadron tensor now reads 

H% = T^(T^-^Bt(D* - Dn)p 3a , P30 ,S aa '(p 2 )SW'(p 2 ) (32) 
z IPs | 

where S aa (p 2 ) — —g aa + P 2 P 2 / m 2 is the standard spin 1 tensor, p 2 = p 3 + Pa, P 3 = Tn 2 -,, pi = m 2 , and p 3 and p± 
are the momenta of the D and the ir, respectively. The relative configuration of the (D,tt)- and (iZ)-planes is shown 
in Fig. HI 

In the rest frame of the D* one has 

P% = (mi,., 6), (33) 

P^ = (E D , |ps| sinr , 0, -|p 3 | cosfl*), 

p% = (E^, -|ps| sinr , 0, |p 3 | cosfl* ), 

|p 3 | = X 1 ' 2 {m%^m 2 D ,ml)/{2m D ,). 

Without loss of generality we set the azimuthal angle \* °f the (D, 7r)-plane to zero. According to Eq. 124|) the rest 
frame polarization vectors of the D* are given by 

e£(±) - -±=(0, ±1, -i, 0), (34) 

4(0) = (0, 0, 0, 

The spin 1 tensor S aa (^2) is then written as 



v , .n 



S aa '( P 2) = -g aa ' + ^f= ^(m)et a '(m). (35) 

Following basically the same trick as in Eq. I|17|l the contraction of the lepton and hadron tensors may be written 
through helicity components as 

L {k)^ H ii v = e »\ m )J v \n)Ll, v ,g mn ,g nn ,e^(m')e v (n')H% (36) 
= L\m,n)g mml g nnl (^(m')el a (r)T; a ) (e^(n')e\ a (s)T^ ^ 

. t 3 Br(Z>* ->£>*■) 
x p 3 e 2 (r) -p 3 e 2 (s) ; 

1 lP3| 



3Br(L>* -> Dn) 



L\t,t)\H^(t)\ 2 -(p 3 el(Q)f 



2|pa| 

£ L k (m,n)H l (m)H^(n)-p 3 e 2 (m)-p 3 ei(r 



m,n— ±,0 



^ L fc (*,n)^(t)ffti( n ) . P3 e 2 (0) .p 3 4W 



n=±,0 



- £ L fc (m ; i)iJ l (m)iJ«W.p 3e2 ( m ).p 3 4(0) 



m=±,0 
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Using these results one obtains the full four-fold angular decay distribution 

dT(B c -> £>*(-> Dir)ll) 



dq 2 d cos 9 d(x/27r) d cos 9 
3 



Br(D* -> Dtt) 



(37) 



Id 



3 . 2 
+ - sm 
4 



-- sm 
4 



cos 2 9) ■ - sin 2 i 
4 



-cos 



2 



drb 1 



drlf 



° ■ i 

cos 2 y • - sm 

4 



+ — sin 20 • cosy • sin 29* 
16 



COS ( 



sin 2 6* 



dq 2 


d<z 2 


r.l 




2 








K 




dr} 2 




d? 2 





dq 2 dq 2 



dVf 



dT 22 



dq 
dTf 
dq 2 



9 

~16 
9 

32 
9 

32 
3 

~4 
3 
2 

3 

~4 



sin 9 ■ cos x • sin 29* ■ 
sin# • sinx ■ sin 29* 
sin 29 ■ sinx ■ sin 29* 
sin 2 6* • sin 2\ • sin 2 (9 



dr\ 2 dr 2 2 



dq 2 dq 2 



dT 12 
aL II | 


dTfA] 


dq 2 


dq 2 )\ 


'dT)\ 


dTf A \ 


v dq 2 


dq 2 J 


(dTf T 


dTf T \ 


\ dq 2 


dq 2 J 



dq 2 8 1 



dT 




dq 2 



IdTf 
4 dq 2 



sin 2 9 ■ sin 2y • sin 2 

16 

Integrating Eq. (|37J) over cos 9* and \ one recovers the two- fo ld ( q 2 , cos 9) distribution of Eq. (|28|l . Note that a similar 
four-fold distribution has also been obtained in Refs.([3,[2|-j2|) using, however, the zero lepton mass approximation. 
If there are sufficient data one can attempt to fit them to the full four-fold decay distribution and thereby extract 
the values of the coefficient functions dTx/dq 2 and, in the case I = t the coefficient functions dTx/dq 2 . Instead of 
considering the full four-fold decay distribution one can analyze single angle distributions by integrating out two of 
the remaining angles as e.g. discussed in Sec. 3. 2 Observables related to single angle distributions will be discussed in 
the next subsection. 



D. Physical observables 



The four-fold distribution Eq. (|37|l allows one to define a number of physical observables which can be measured 
experimentally. An asymmetry parameter ag» is defined from the angular distribution W(cos 2 9*) = 1 + ag* cos 2 9* . 
Integrating Eq. Ij37(l over cos 9 and x one finds 
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-{U 11 + U 22 ) + 2{L 11 + L 22 ) - U 11 + 2U 22 + 2(L n - 2L 22 ) + 2S 22 

Ug* = = (38) 

U 11 + U 22 + U 11 - 2U 22 

By integrating over cos 9* and \ one can define two asymmetry parameters a e and ag according to the angular 
distribution W (cos 2 9) — 1 + a' e cos 9 + ag cos 2 9. One has 



Ijn + v ii + 2 (£H + L 22 ) + 2(U n - U 22 - 2L 22 + (2/3)S 22 ) 



(39) 



[/ll + _ 2 ( L 11 + L 22) _ 2(£/ ll + ^22) + 4(L 11 + L 22) 

" e ~ U 11 + U 22 + 2(L n + L 22 ) + 2(U 11 - U 22 - 2L 22 + (2/3)S 22 ) ' 

An azimuthal asymmetry parameter (3 can be defined from the ^-distribution 
W(x) = 1 + /3cos2x- The azimuthal ^-distribution is obtained by integrating over cos 9 and cos 9*. One has 

3 = -(T ll +T 22 ) + 2(f 11 +f 22 ) 

~ U 11 + U 22 + L 11 + L 22 + U 11 - 2U 22 + L 11 - 2L 22 + S 22 ' 

A second strategy is to define suitable asymmetry ratios that project out the partial rates from Eq. (|37fl . Let us 
consider the following four asymmetry ratios which project out the contributions of the parity conserving partial rates 
Tt and Tj and the parity violating partial decay rates Tp and Ta- One has 

r A _ dT(x) - dT( X + it/2) + dT( X + tt) - dT( X + 3tt/2) 

T - T dr(x) + dr( x + 7r/2) + dr(x + 7r) + dr( x + 3^/2) 1 J 

where — 7r/4 < x < 7r/4. 

T A : A I= Ni/Dj , (43) 



Ni = dT(e,e*, X )-dT(8,d*,x + n) 

- dv(e,-n:-e*,x) + dT(e,Tr-e*,x + ^) 

- dr{ir-9,9*,x)+dT{Tr-9,9*,x + n) 

+ dT(Tr - 0, tt - 9*, x) - dT{ir -e,TT-9*,x + 7r), 

7T 7T 7T 7T 

< 6»* < - , _<6»<7T, <X<~- 

- ~2' 2~-' 2~ A ~2 

The denominator Dj is given by the same expression with plus signs everywhere. 

dT(6) — <£T(tt — 6) tt 

. dT(g* , X ) ~ dr(9*, X + ir) - dT(n - 9*, X ) +dT(n - 9*, X + tt) 

A ' A dT(9*,x)+dT(9*,x + K)+dT(Tr-9*,x)+dT(Tr-9* 1 x + Tr) 1 

TT TT TT , 

0<^<-, __< x <_. (45) 

We have used a notation where the angles that do not appear in the arguments of the differential rates dT have been 
integrated out over their physical ranges (0 < 9 (9*) < tt, < x < Itt). Integrating over the remaining variables 
(numerator and denominator separately!) we finally obtain 
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cAt . 2 -(i/2)(r" + r") + r" + r" ) 
7r r 

2 (l/2)(I 11 +/ 22 )-7 11 -/ 22 

7T 1 

3?; P 12 

<^ B > = — — , (48) 
3vA 12 +A 21 

<A A > = — . (49) 

E. Longitudinal polarization of the lepton in B — *• A7/-decay. 

Our aim is to study lepton polarization effects in B — > Kll-decsy. The longitudinal polarization of the final lepton 
l~{k2) is defined by 

p(0 = rfr( S2 )/dg 2 - dT{-s 2 )/dq 2 

dT{s 2 )/dq 2 +dT{~s 2 )/dq 2 [ ' 

where the longitudinal component of the rest frame polarization vector of the l~ (k 2 ) is given by 

Contrary to previous studies of the longitudinal polarization of the lepton, where the longitudinal polarization was 
studied in the B rest frame, we will calculate this quantity in the ll-CM frame. The longitudinal polarization vector 
s 2 is boosted to the moving frame by a Lorentz transformation. One obtains 

f|k 2 | E 2 k 2 \ 

s 2 ,cm = , — Tj — r • (52) 

V A 1 M l k 2|/ 

The quantity dT(s 2 ) / dq 2 in Eq. (|50|l may be obtained from Eq. I|14|) by the replacement 

(fc+V)^ ^(1 + 75 &)(#2+M). 

Integrating the numerator and denominator in Eq. (|50|l over cos 9, one finds 

m v (dU 12 + dL 12 ) , 

p = — df (53) 

where we have adopted a short hand notation U :— Tjj, U := Tjj, etc.. Because of CP-invariance and because 
the longitudinal polarization is a pscudoscalar quantity, the longitudinal polarization of the antilepton is equal and 
opposite to the longitudinal polarization of the lepton, i.e. P^(l + ) = — pW(i~). 



IV. MODEL FORM FACTORS 



We will employ the relativistic constituent quark model jTsL IT^ | to calculate the form factors relevant to the decays 
B — > Kll and B c — > D(D*)ll. This model is based on an effective interaction Lagrangian which describes the coupling 
between hadrons and their constituent quarks. 

For example, the coupling of the meson H to its constituent quarks qi and q 2 is given by the Lagrangian 

C mt {x) = g H H{x) Jdx 1 Jdx 2 F H (x,x 1 ,x 2 )q(x 1 )r H X H q(x 2 ) . (54) 
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Here, A# and Th are Gell-Mann and Dirac matrices which entail the flavor and spin quantum numbers of the meson 
field H(x). The function Fh is related to the scalar part of the Bethe-Salpeter amplitude and characterizes the finite 
size of the meson. The function must be invariant under the translation Fs(x + a, x\ +a, X2 + a) = Fh(x,Xx,x 2 ). 

In our previous papers we have used the so-called impulse approximation for the evaluation of the Fcynman diagrams. 
In the impulse approximation one omits a possible dependence of the vertex functions on external momenta. The 
impulse approximation therefore entails a certain dependence on how loop momenta are routed through the diagram 
at hand. This problem no longer exists in the present full treatment where the impulse approximation is no longer 
used. In the present calculation we will use a particular form of the vertex function given by 



F H (x,Xx,x 2 ) 



m\Xi + m 2 x 2 
mi + m 2 



(55) 



where m\ and m 2 are the constituent quark masses. The vertex function Fh evidently satisfies the above translational 
invariance condition. As mentioned before we no longer use the impulse approximation in the present calculation. 

The coupling consta nts q h in Eq. (|54|) are determined by the so called compositeness condition proposed in [2(j 
and extensively used in |21j . The compositeness condition means that the renormalization constant of the meson field 
is set equal to zero 



4tt 2 



U' H (m H )=0 



(56) 



where II ^ is the derivative of the meson mass operator. For the pseudoscalar and vector mesons treated in this paper 
one has 



2p 2 dp a J Air 2 i 



x tr 

1 

3 

x tr 



p^p v 



where tuy = mj/(mi + mj), <&#(— A; 2 ) is the Fourier-transform of the correlation function ((xi — x 2 ) 2 ) and Si(ji) 
is the quark propagator. The lcptonic decay constant fp is calculated from 

3gp f d A k 



47T 2 / 47T 2 



fpp" 



(57) 



The transition form factors P(px) — * P(p 2 ), V{p 2 ) can be calculated from the Feynman integral corresponding to the 
diagram of Fig. [5] 



k T (px,p 2 ) 



3gpgp>(V) f d A k ~ 



4tt 2 



47T 2 



■$p(-(fc + wxapx) ) $p'(v)(-{k + w 23P2 y) 



x tr 



S 2 (ft+ h)^Sxm+ ^ih 5 S 3 (]i)T out 



(58) 



where = 7 M , 7 M 7 , ia^"q v , or ia^q^ 5 and T pi y — 7 , 7w e 2- 
We use the local quark propagators 



mi— ft ' 



(59) 



where mj is the constituent quark mass. We do not introduce a new notation for constituent quark masses in order to 
distinguish them from the current quark masses used in the effective Hamiltonian and Wilson coefficients as described 
in Sec. II because it should always be clear from the context which set of masses is being referred to. As discussed in 
UM E3 > w e assume that 



mH < mx + m 2 



(60) 
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in order to avoid the appearance of imaginary parts in the physical amplitudes. 

The fit values for the constituent quark masses are taken from our papers [TsUli^ and are given in Eq. (|61[) . 

m u m s m c m b 

(61) 



0.235 0.333 1.67 5.06 GeV 

It is readily seen that the constraint Eq. (|60|l holds true for the low-lying flavored pseudoscalar mesons but is no 
longer true for the vector mesons. In the case of the heavy mesons D* and B* we will employ identical masses for 
the vector mesons and the pseudoscalar mesons for the calculation of matrix elements in Eqs. (|56|) . l|57l) and (|58() . It 
is a quite reliable approximation because of (tbd* — mo)/m£i ~ 7% and (mj. — mflj/m^ ~ 1%. 

We employ a Gaussian for the vertex function {k 2 E / A 2 H ) — exp(— k 2 E / A 2 H ) where fee is the Euclidean momentum 
and determine the size parameters Ah by a fit to the experimental data, when available, or to lattice simulations 
for the leptonic decay constants. The quality of the fit can be seen from Table IIIII The branching ratios of the 
semileptonic decays are shown in Table Hvl The numerical values for Ah are A ff = 1 GeV, Ak = 1.6 GeV, Ad = 2 
GeV and A E = 2.25 GeV for all K, D and B partners, respectively. 

We are now in a position to present our results for the B — > K form factors. We have used the technique outlined 
in our previous papers [l8L Il9| for the numerical evaluation of the Feynman integrals in Eq. I|58l) . The results of our 
numerical calculations are well represented by the parametrization 

(62, 

1 — as + bs 

Using such a parametrization facilitates further integrations. The values of ^(QV, a and b are listed in Tables IVIVII 
We plot our form factors in Fig. [3] and compare them with those used in paper |8( in Fig. The functional behavior 
of the curves is similar to each other. 

At the end of this section we would like to discuss the impulse approximation used in our previous papers [Tsl Il9| . 
It was simply assumed that the vertex functions depend only on the loop momentum flowing through the vertex. The 
explicit translational invariant vertex function in Eq. I|55[) allows one to check the reliability of this approximation. 
We found that the results obtained with and without the impulse approximation are rather close to each other except 
for the heavy-to-light form factors. We consider the B — > 7r-transition as an example to illustrate this point. The 
calculated values of the Ff 7r (q 2 ) form factor at q 2 = are 



0.27 exact 



0.48 impulse approximation 



One can see that the value of the form factor at q 2 = calculated without the impulse approximation is considerably 
smaller than when calculated with the impulse approximation. Its value is close to the value of QCD SR estimates, 
see, for example, [2|j: 



Ff w (0) = 



0.25 asymptotic distribution 
0.30 QCD SR distribution 



V. NUMERICAL RESULTS 



In this section we collect and discuss our numerical results. We plot the normalized differential distributions 
T^ \ dT/ds with r to t = 1/tb (tb = (r B o + t b +)/2 — 1.60 ps) and s — q 2 /m 2 B in Figs. 15161 for the decay B — > Kll. 
We have also included the B — ► Kvv modes. Their differential rates are calculated according to 



dT (B+ -> Kuu) = G 2 F m%\X t \ 2 a 2 \D u {x t )\ 2 
ds 2 8 7T 5 sin 4 6 W 



,2 2 



F+\4>1- (63) 



The functions D v (xt) and <f>H are given in Appendix B. The behavior of the normalized differential distributions is 
shown in Fig. [7| 

We list our numerical results for the branching ratios in Table Ivnl When comparing the values of the branching 
ratios with those obtained in [j| one finds that they almost agree with each other. 
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Finally, we plot the dependence of the normalized differential distributions s = q 2 /m 2 B in Figs. 1101121 for the 
decay B c — > D(D*)Tl(pv) . In the numerical analysis we use the input parameters: itlb c — 6.4 GeV, tb c — 0.46 
ps and IV^Vtbl = 0.008. The B c — > /^(D^-transition form factors are plotted in Figs. 15151 and the normalized 
differential distributions for B c D(D*)fj,+^-, B c -> D{D*)t+t- and B c D(D*)vv are shown in Figs. I10I11IL21 
respectively. The results for the branching ratios are also given in Table IVIII They are to be compared with the 
results of calculations performed in [Tt| where the light front and constituent quark models were employed. 
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Appendix A: Wilson Coefficients 



In this paper we use the Wilson-coefficients d calculated in the naive dimensional regularization (NDR) scheme in 
the leading logarithmic approximation Q: 



CjQi) = J2 k ^ C7' = l,».6) 

i=i 

8 8 
Clip) = V^C 7 (M W ) + - -t^J C s (M w )+J2hiV ai , 



i=i 



with 



a(M w ) 
a(jjt) 



C 7 (M W ) = -~A(x t ), C 8 (M W ) = ~F(x t ), 



where xt = mf/M^ and A(x) and F{x) are defined below. The numbers a%, fcjj and hi are given in Table IVTTT1 
The coefficient of Qio is given by 



C W (M W ) 



Y(x t ) 
sin 2 6» 



with Y(x) given below. Since Qio is not subject to renormalization under QCD, its coefficient does not depend on 
/i ps 0(rrib). The only renormalization scale dependence enters through the definition of the top quark mass. 
Finally, including leading as well as next-to-leading logarithms, one finds 

C 9 (ji) = P + J[ X ^ - 4Z{x t ) + P E E(x t ) 



with 



sin 6 1 



P ° = fir \ (-0- 1875 + Y,PM ai+1 ) + 1-2468 + ]T V a > [n + s iV ] 
a s (M w ) ^ ^ 



i=i 



Pe = O.U05 + J2^V az+1 



Here 



Y(x) = C(x) - B(x), Z{x) = C(x) + -D(x). 



A(x 
B(x 

C(x 

D(x 

E(x 

F(x 



a; (8a; 2 + 5a; - 7) x 2 (2-3x) 
12(x- l) 3 + 2(x-l) 4 

In a;, 



hi . 



4(1 -x) A(x-1) 2 
x(x — 6) x(3x + 2) 



In a;, 



(x-1) 8(a;-l) 2 
-19x 3 + 25x 2 x 2 {5x 2 - 2x - 6) 



18(*-1)« lna; -9 lnX ' 



36(x-l) 3 
x(X8 - 11a; - x 2 ) a; 2 (15 - 16a; + 4a; 2 ) 



12(1 -x) 3 
a;(a; 2 — 5a; — 2) 



A(x- l) 3 2(ar- l) 4 
The coefficients Pi, Sj, and are given in Table IVIIII 



6(1 -x) 4 
In a;. 



In a; In a;, 
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Appendix B: The functions in Eqs. (|31|) and (1631) . 

We list here a set of the functions appearing in Eqs. EU and (gSJ from 8] and [l7|. 



D v (x) 



l-r H j -2sl l + rj* i i-.s" 
x (2 + x 3x — 6 



35—1 (x — 1) 
2 



Inx 



ai = (1 - \/th) 2 

(i - V^) 2 



(i + V^7) 2 



/3i = 



If 1-8 



A 



V 



V 



- 1 



(i + 

1 ^ima^) 







2 


1r H {l 


+v 





Q 2 



2 \ r H J1 + V^" 
ff , 2m b CfF T 2 

2 



\C W F+ 



P2 = \C W \ 

"3 = (1 - Vth) 2 



1 + V^" 

1 + r„ - |F + | 2 + (l - r H ^jRc(F + Fl) + i S |F_ 



ff 2m b Cf(l + y?F)ao 
G 9 H 



ft 



(1 + v^i 7 ) 2 
s 

(1 + V^ 7 ) 2 



cfv + 



2m b Cf{l + ^)g 



2 




2-i 


+ 


CioAq 






2-i 





c w v 



s 2m b Cf (l + VFF)ao 

Og A H 



2 




2-i 


+ 


CioA) 





to 




2-i 


+ 







2rh b Cf (1 + VrF)a 



ifi^-l)i^Re 

2 \ r H y 1 + v^F 



CfA 



5 = 



CfA + 

ICiol 2 
2(1 + V^F) 2 

(j>HS 



2m b Cf ff (l + V ?F) a " 
s 

|Cio| 2 Re(A4) 



-20 ff |y| 2 -3(l-r ff ) 2 |A | 5 



2(l + r H )-s 



I Ah 



4r 



^ | A _|2 + ;^___iIZ R e( _ A A\ - A Al_ + A+Al ) \ . 



where rh b = m b /m,B- 
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TABLE I: Central values of the input parameters and the corresponding values of the Wilson coefficients used in the numerical 
calculations. 



mw 


80.41 GeV 


Ci 


-0.248 


mz 


91.1867 GeV 


C 2 


1.107 


sin 2 6w 


0.2233 


c 3 


0.011 


m c 


1.4 GeV 




-0.026 


rot 


173.8 GeV 


C 5 


0.007 


?7lb,pole 


4.8 GeV 


a 


-0.031 




?^b,pole 


cf 


-0.313 




0.220 GeV 


c 9 


4.344 


a" 1 


129 


Cio 


-4.669 


a 8 {m z ) 


0.119 


Co 


0.362 


\v t \v tb \ 


0.0385 






\v? d v tb \ 


0.008 






\v t \v tb \/\v cb \ 


1 







TABLE II: Bilinear combinations of the helicity amplitudes that enter in the four- fold decay distribution Eq. ( 1371 '). 



Definition 


Property 


Title 


H}i = Re 
H% = Im 




+ Re| 
) +Im 


{hIhP^ 


) 


rrij _ rrji 
n V ~ n U 

rrij _ _ zrji 
n iv — n lV 


Unpolarized- 
transverse 


H*J = Re 1 
H% = Im 




-Re| 
-Im 


HtHp) 
(hLhP 


) 


H% = H P ' 

ttH ttH 

H/ P = -H 3 IP 


Parity-odd 


Bi} = Re I 
H% = Im 


'H\Hp) 






Transverse- 
interference 


H$ = Re ( 
H\i = Im 






HI = H£ 
zrij — _ui i 

n IL — n lh 


Longitudinal 


Hg = 3 Re (HlHj j 


) 
) 


Hi j = hi; 

Tjij _ _ rrji 

n is — n IS 


Scalar 


fl& = Re 
H% L = In 




) 




Scalar- 
Longitudinal- 

interference 


H? = | [] 
H\{ = \ \ 


[m ( 


H\H 
W+H 


V) 


+ Re ( 
+ Im( 


HlHp)\ 






transverse- 

longitudinal- 

Interference 


HI = | [] 
fljji = § \ 


}e ( 
Im 




-Re( 
-Im 


[hLhP) 


[ 




parity- 
Asymmetric 


Tjij _ I 
"ST — 2 


Re 
|"lm 




+ Re 
) +Im 


[lit Hp) 

(hlhV 


i 




Scalar- 

Transverse- 

interference 


-"SA — 2 

H ij — i 
■"ISA — 2 


Re 
Tim 


(h'+hP 


- 

)- 


Re 
-Im 


(hihP) 

(Hi Hp 


>i 




Scalar- 

Asymmetric- 

interference 



19 



TABLE III: Leptonic decay constants ju (MeV) used in the least-square fit. The values are taken either from PDG [30| or 
from the Lattice |31| : quenched (upper line) and unquenched (lower line). 



Meson 


This model 


Expt/Lattice 


7T+ 


131 


130.7 ±0.1 ±0.36 


K+ 


161 


159.8 ± 1.4 ±0.44 


D+ 


211 


203± 14 






226± 15 


Df 


222 


230± 14 






250± 30 


B+ 


180 


173± 23 






198± 30 


B° s 


196 


200± 20 






230± 30 


B+ 


398 





TABLE IV: Semileptonic decay branching ratios. 



Meson 


This model 


Expt. 


n + _> ^l+v 


1.03 • 10 -8 


(1.025 ± 0.034) ■ 10~ 8 


K + 


4.62 ■ 10~ 2 


(4.82 ±0.06) ■ 10" 2 


B + ^D°l+u 


2.40 • 10 -2 


(2.15 ±0.22) ■ 10~ 2 


B + ^D*°l+v 


5.60 ■ 10~ 2 


(5.3 ±0.8) ■ 10" 2 


B+ -» D°l + u 


2.05 • 10~ 5 




B + ^D*°l+v 


3.60 ■ 10" 5 





TABLE V: Parameter values for the approximated form factors F(s) = F(0)/(1 — as + bs 2 ) (s = q 2 /m 2 B ) in the decays B — * Kll. 





F+ F- F T 


a 
b 


0.357 -0.275 0.337 
1.011 1.050 1.031 
0.042 0.067 0.051 



TABLE VI: Parameter values for the approximated form factors F(s) — F(0)/(1 — as ± 6s 2 ) (s = q 2 /m%) in the decays 
B c D(D*)Tl. 





F+ F- F T 


A A+ A- V 


fflo a+ g 


no) 

a 
b 


0.186 -0.190 0.275 
2.48 2.44 2.40 
1.62 1.54 1.49 


0.279 0.156 -0.321 0.290 
1.30 2.16 2.41 2.40 
0.149 1.15 1.51 1.49 


0.178 0.178 0.179 
1.21 2.14 2.51 
0.125 1.14 1.67 
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TABLE VII: Decay branching ratios without (with) long distance contributions. 



Ref. 


Br(B -> Ky, + n~) 


Br(B -> Kt + t~) 


Br(B — > K vv) 


[8] 


0.57- 10" 6 


1.3- 10" 7 




\£ 


(0.35 ±0.12) ■ 10" 6 






Ul] 


0.44 • 10 -6 


1.0- 10~ 7 


5.6- 10~ 6 


[32] 


0.5 • 10~ 6 


1.3- 10" 7 




our 


0.55 (0.51) ■ 10" 6 


1.01 (0.87) ■ 10" 7 


4.19 ■ 10~ 6 





our 


[12] 


Br(B c -> D d /i+^i") 

Br(B c -> I^/i+zO 
Br (B c -> £> s (j+nl 
Br(B c -> D*fj, + n-) 


0.44(0.38) ■ 10 -8 
0.71 (0.58) ■ 10" 8 
0.97(0.86) ■ 10" 7 
1.76(1.41) ■ 10" 7 


0.41 (0.33) ■ 10" 8 
1.01 (0.78) ■ 10~ 8 
1.36(1.12) • 10" 7 
4.09 (3.14) ■ 10" 7 


Br(B c -> D d r+r-) 
Br(S c -^D* d r+T-) 
Br (B c -> £> s t+t~) 
Br (B c -> D:r + r-) 


0.11(0.09) ■ 10" 8 
0.11(0.08) ■ 10~ 8 
0.22 (0.18) ■ 10~ 7 
0.22 (0.15) • 10~ 7 


0.13(0.11) ■ 10 -8 
0.18(0.13) ■ 10~ 8 
0.34 (0.27) ■ 10~ 7 
0.51 (0.34) • 10" 7 


Br(B c -> L> d Pi/) 
Br(B c -> D* d vv) 
Br(B c -> D s vv) 
Br(B c -> D* a vv) 


3.28 • 10" 8 
5.78 • 10" 8 
0.73 ■ 10~ 6 
1.42 • 10" 6 





TABLE VIII: Values of parameters in the formulae for the Wilson coefficients. 





1 


2 


3 


4 


5 


6 


7 


8 


ai 
ku 

k 3i 

kii 

k 5 i 
km 


14 
23 










16 
23 










6 

23 
1 
2 
1 
2 

1 

14 
1 

14 





12 
23 
1 
2 
1 
2 
1 
6 
1 
6 





0.4086 



0.0510 
0.0984 
-0.0397 
0.0335 


-0.4230 



-0.1403 
0.1214 
0.0117 
0.0239 


-0.8994 



-0.0113 
0.0156 
-0.0025 
-0.0462 


0.1456 



0.0054 
0.0026 
0.0304 
-0.0112 


hi 
h 


2.2996 
0.8623 


-1.0880 



3 
7 




1 

14 




-0.6494 
-0.9135 


-0.0380 
0.0873 


-0.0185 
-0.0571 


-0.0057 
0.0209 


Pi 
r% 
s t 
1i 














80 
203 

0.8966 
-0.2009 



8 

33 

-0.1960 
-0.3579 



0.0433 
-0.2011 
0.0490 
0.0318 


0.1384 
0.1328 
-0.3616 
0.0918 


0.1648 
-0.0292 
-0.3554 
-0.2700 


-0.0073 
-0.1858 
0.0072 
0.0059 




FIG. 2: Diagram describing the form factors of the decay B(B C ) — > K(D, D*)ll 
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0.2 0.4 0.6 0.8 

s 

FIG. 7: Normalized differential distributions 10 s F^dT/da for B -> Kvv. 




0.1 0.2 0.3 0.4 0.5 



s 

FIG. 8: Form factors for the B c -» D d transition: (a) F + , (b) F-, (c) Ft. 
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FIG. 9: Form factors for the B c — > transition. Upper panel: (a) V, (b)Ao, (c) A+, (d)A-. Lower panel: (a)g, (6) do, 
(c)o+. 
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